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ABSTRACT 

We use cosmological hydrodynamic simulations to consistently compare the assembly of 
dwarf galaxies in both A dominated, Cold (CDM) and Self-Interacting (SIDM) dark mat¬ 
ter models. The SIDM model adopts a constant cross section of 2 cw? jg, a relatively large 
value to maximize its effects. These are the first SIDM simulations that are combined with a 
description of stellar feedback that naturally drives potential fiuctuations able to create dark 
matter cores. Remarkably, SIDM fails to significantly lower the central dark matter density 
within the central 500pc at halo peak velocities ^max < 30 km s“^. This is due to the fact that 
the central regions of very low-mass field halos have relatively low central velocity dispersion 
and densities, leading to time scales for SIDM collisions greater than a Hubble time. CDM 
halos with ^max < 30 km s“^ have inefficient star formation, and hence weak supernova feed¬ 
back. At a fixed 2 cm? jg SIDM cross section, the DM content of very low mass CDM and 
SIDM halos differs by no more than a factor of two within 100-200pc. At larger halo masses 
10^^ Mq), the introduction of baryonic processes creates field dwarf galaxies with dark 
matter cores and central DM+baryon distributions that are effectively indistinguishable be¬ 
tween CDM and SIDM. Both models are in broad agreement with observed Local Group field 
galaxies across the range of masses explored. To significantly differentiate SIDM from CDM 
at the scale of faint dwarf galaxies, a velocity dependent cross section that rapidly increases 
to values larger than 2 cw? jg for halos with ymax < 25-30 km s“^needs to be introduced. 
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1 INTRODUCTION 

In this paper we use high resolution SPH + N-Body cosmological 
simulations to focus on two main questions: Can we identify a clear 
difference between Cold Dark Matter and Self-Interacting Dark 
Matter (CDM and SIDM respectively) predictions for the structural 
and observable properties of dwarf galaxies? How do the properties 
of CDM and SIDM halos differ at halo masses below 10^° M©, 
where star formation (SF) becomes very inefficient and the effect 
of the underlying DM component should dominate? 

SIDM was originally introduced as a solution for the so called 
“core-cusp” problem, the excess of dark matter predicted by the 
CDM model at the center of field and satellite dwarf galaxies com¬ 
pared to observations ([Moore et al.|1999bl|Oh et al.|2008|| Walker | 
|& Penarrubia|201l|[Adams et al.|2014| )7hi the SIDM model dark 

matter particle collisions isotropize the cores of galaxies and trans¬ 
fer mass outward from the dense central regions of DM halos over 
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cosmic time scales. This process creates large cores, more spherical 
halos, and a signature flat radial profile of the DM velocity disper¬ 
sion ( [Spergel & Steinhardt[[2000|[Burkert[[2000[ ). A SIDM model 
with fixed cross section and elastic collisions represents the sim¬ 
plest model in a large class of plausible “dark sector” DM models. 
However, significantly more complex interactions are possible, in¬ 
cluding Yukawa potentials ( [Feng et al.[2009[[20r0|[Loeb & Weiner 


2011 1, and cooling or atomic dark matter ([Cyr-Racine & Sigurdson 


2013 [Schutz & Slaty er|2015[ [Buckley et al.|2014| |. The dynamics 

of SIDM was first implemented in cosmological simulations as a 
fluid ( [Moore et al.[2000[| and then as elastic collisions between par- 


ticles ([Spergel & Steinhardt|2000[[Burkert|2000 pave et al.|2001 


Colin et al.[2002[[Strigari et al.[2007[[Martinez et al.[2009[ Koda &[ 
Shapiro|201 l[[Vogelsberger et al.|2012[l~ 


Numerical studies combined with observations (the elliptic- 
ity of galaxy clusters, the abundance of and the core size of held 
dwarf galaxies) have constraine d the cross section of SIDM to be 
of the order of ~ 0.1 -1 cm ^/g |Loeb & Weiner 2011 Rocha et al. 


[2013| [Peter et aT][2013| [Zavala et al.[[2013[ l. However, a variable 


© 2012 RAS 































2 Bastidas-Fry et al. 


SIDM cross section (a decreasing function of the DM particles’ 
relative velocity) can preserve the ellipticity and density of galaxy 
clusters while producing cores in dwarf galaxies ( [Peter et al.|2013| 
but see [Newman et al.j ( [2013[ ) for evidence for DM cores in clus¬ 
ters). While weakening existing constraints, a variable SIDM cross 
section is well motivated by “hidden sector” particle physics mod¬ 
els ( [Tulin et al.|2013[l. Further c rucial constraints on large SIDM 
cross sections ( [Elbert et al.[20I4[ ) come from the necessity of form¬ 
ing cores in very faint galaxies without evaporating the satellites of 
MW-like halos and galaxy sized halos in clusters as they interact 
with a dense DM environment ( [Gnedin & Ostriker|2001[ [Vogels- [ 
[berger et al.[ 2012 [ ). 

Unfortunately, a large fraction of the astrophysically driven 
support for non-standard DM models has come so far from simpli¬ 
fied simulations that lack the complexities of “baryon physics” and 
follow only the the assembly of the DM component (see the review 
by [Brooks[2014[ ). The necessity to couple a DM model with baryon 
physics comes from the existence of “bulgeless galaxies”, a prob¬ 
lem that requires the removal of low angular momentum gas from 
galaxies (Binney et al.[200l|[van den Bosch et al.[2001| Governato| 


[et aL|2010| l through feedback processes ( Brook ct aL|2011[ l. Fur- 

ther motivation for including baryon physics includes the necessity 
to que nch SF in galaxy satellites (|Klypin et aL|19^ [Moore et~^ 

[1999a[ [D’Onghia & Burkert|2003| , especially if the SIDM power 

spectrum is similar to CDM at small scales and subhalos survive 
evaporation (as found by [Zavala et al.[2013] l. 

Crucially, analytical and numerical work have shown that 
feedback lowers the central DM density in galaxies, creating gas 
outfiows and repeated fiuctuations in the gravitational potential 


([Mashchenko et al. [2008 [[Del Popolo[2009 

Governato et al.[2010| 

Pontzen & Govemato|2012| Di Cintio et a] 

l.|2014|[Teyssier et al.[ 

2013|[Velliscig et al.[2014[l. This results in irreversible energy trans- 


fer to the DM (see also review by [Pontzen & Governato[[2014| l. 
These outflows, generated by a bursty SF, have strong observational 
support ( van der Wei et al.[ 20 TT [Martin et al.[ 2012 |[Lundgren et al.[ 


[2012[ [Dominguez et al.[2014[[Geach[2015[ l. Bursty SFH have also 

been robustly associated with the b uild up of the stellar content of 


galaxies in the 10^-10^° Mprange ( [Tolstoy et al. 2009 Kauffmann 


[2014[ ). The ability of feedback to dynamically heat the DM can po¬ 
tentially remove the need for SIDM at galactic scales. As a result 
it is fundamental to identify the unique differences between CDM 
and SIDM when both models include an explicit treatment of the 
physics of galaxy formation. 

This study is able to follow the evolution of field SIDM halos 
at masses below ~10^° Mq (VmaaQ< 40 km s“^). This is an im¬ 
portant regime as ( 1 ) observational data are becoming robust (e.g., 
[Papastergis et al.[[2015[ ), and (2) feedback processes become less 
efficient with declining Yma x ([Governato et al. 12012) [Penarrubiaj 
[et al.[[201^ [Pi Cintio et al.[[^14[ ). The most recent simulations 

( [Madau et al.[[2014] [Governato et al.[[T015| [Onorbe et al.[[2015T ) 

show core formation when only 10 ^'^“® M 0 of stars have formed. 
Also, it is often naively assumed that SIDM will form DM cores 
at all dwarf masses. However, at these small halo masses the trend 
of later assembly epoch of the first prog enitor ([Li et al.[2008[l, and 
consequentially lower central densities ( [Avila-Reese et al.[[2005| l, 
will affect the timescale and extent of SIDM core formation in field 
galaxies beyond simple scaling calculations. Furthermore, the rela¬ 
tive velocity of a halo and the surrounding DM background is cru- 


^ where Nmax is defined as the peak of the rotation curve, with Nrot = 

(GMlrYl^. 


cial in determining the collision rate, boosting it in small halos (that 
have intrinsic low a). Our subset of DM-only simulations is also 
one of the first to compare very small field and satellite galaxies at 
similar mass and spatial resolution. The central density of field and 
satellite dwarfs (that move fast through a dense DM environment) 
may then be significantly different, but it has not been properly 
compared yet. 

In this work we adopt a constant velocity cross section of 2 
cm^/g in the SIDM runs (a relatively large value to maximize its ef¬ 
fects) and a common description of SF and feedback in both CDM 
and SIDM models. The SF and feedback prescriptions have been 
shown to form CDM galaxies with SF efficiency, photometric, and 


kinematic properties close to those of real ones ([Oh et al.[ 

12011 

Munshi et al.[[2013| [Christensen et al.[[2014a| [Shen et al.| 

2014 

Brooks & Zolotov|2014|). These simulations are the first high res- 


olution simulations to compare the assembly of dwarf galaxies in 
CDM vs SIDM cosmologies including a description of SF and su¬ 
pernova feedback that creates realistic galaxies while creating DM 
cores through “DM dynamical heating.” Other recent numerical 
work ( [Vogelsberger et al.[2014[ | focused on small galaxies in SIDM 
and explored the role played by a variable cross section, however, 
the feedback recipe implemented in their work does not generate 
DM cores. 

In §2 we describe the simulation and the code used. In §3 we 
present the results, which are then discussed in §4. 


2 METHODOLOGY 


2.1 ChaNGA, star formation and bursty feedback with 
outflows 

The simulations were run in a full cosmological context and to 
a redshift of zero using the N-body Treecode + Smoothed Par- 
ticle Hydrodynanfics (SPH) code Cha NGa ( [Jetley et~^[2008[ 

ChaNGa includes stan- 


Quinn et al.|2013| Menon et al.|2014p 


dard physics modules previously used in Gasoline ( [Wadsley et al.[ 
[2004[ [2008[ [Stinson et al.[[2012} including a treatment of metal 
line cooling, self shielding, cosmic UV background, star formation, 
“blastwave” SN feedback and thermal feedback from young stars 
( [Stinson et al.[2006j[2012[). The SPH implementation includes ther- 
mal diffusion ( [Shen et al.[201^ and eliminates artificial gas surface 
tension through the use of a geometric mean density in the SPH 
force expression ( Ritchie & Thomas|2001[[Governato et al.|2015[ 


[Menon et al.[[2014[ t. This update better simulates shearing flows 

with Kelvin-Helmholtz instabilities. For consistency with our pre¬ 
vious work comparing CDM and WDM scenarios, we adopted the 
same feedback and SF pa rameters as inlGovernato et al.[ ( [2015] ). A 
Kroupa IMF is assumed ( [Kroupa|200lj ) and the density threshold 
for SF is set at 100 amu/cm^. Limiting star formation to dense gas 
regions is a realistic approach and concentrates feedback energy 
([Brook et al. [[2012| [Christensen et al.[[2014b| [Agertz & Kravtsov[ 

|2014> . 100% of SN energy is coupled to the surrounding gas. 


^ ChaNGa is part of the AGORA group, a research collaboration with 
the goal of compare the implementation of hydrodynamics in cosmolog¬ 
ical codes <Kim et al.|[20T4t . ChaNGa is available here: http://www- 
hpcc.astro.washington.edu/tools/changa.html 
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Name 

Physics 

DM/gas particle mass Mq 

Force Softening (pc) 

halo mass range (Mq ) 

h516CDM 

SF and DM-only 

1.6 104/3.3 103 

86 

10® - 5 X 1040 

h516SIDM 

SF and DM-only 

1.6 104/3.3 10® 

86 

10® - 5 X lOio 

h2003CDM 

SF and DM-only 

0.67 104/1.4 103 

64 

4 X 10® - 1 X IQi® 

h2003SIDM 

SF and DM-only 

0.67 104/1.4 103 

64 

4 X 10®-1 X 1040 

40Thieves-CDM 

DM-only 

0.81 104 

64 

4 X 10®-2.7 X IQio 

40Thieves-SIDM 

DM-only 

0.81 104 

64 

4 X 10®-2.7 X IQio 

40Thieves-SIDM.hr 

DM-only 

0.24 104 

64 

1.25 X 10® -2.7 X IQio 

hl48CDM 

DM-only 

1.93 104 

86 

10®-2 X 1012 

hl48SIDM 

DM-only 

1.93 104 

86 

10®-2 X 1012 


Table 1. Simulations and physics parameters explored in our simulations. All SIDM runs adopt (jsidm = 2[cm^ g“^]- The mass range 
shows halos with at least 50,000 DM particles within their virial radius. The main halos in each volume are studied with several million particles each. The 
DM-only subsample is one of the first to simulate a sample of field and satellite halos at similar mass and spatial resolution. This approach allowed us to study 
the effect of the environment on the DM distribution of the halos in the two populations. 


2.2 SIDM Implementation and Analytical Expectations 

The SIDM implementation closely follows the standard Monte 
Carlo approach described in previous works Pave et aL|200T|[Wl 
Igelsberger et al.|20T^ , and here is only briefly summarized, while 
tests are presented in Appendix A. SIDM interactions are modeled 
under the assumption that each simulated DM particle represents a 
patch of DM phase-space density and that the probability of colli¬ 
sions is derived from the collision term in the Boltzmann equation. 
Collisions are then elastic and explicitly conserve energy and mo¬ 
mentum. For a detailed discussion see ( [Rocha et al.|2013| l and also 
(jYoshida et al.|2000||D’Onghia & Burkert|2003[|Kaplinghat et al.j 

|2014j t. The number of interactions N that occur in a region with 
local DM density p in a time At is 

N = pvadm^t (1) 

In practice to simulate interactions for discrete N-body DM 
particles we use equation one as the implied probability for a par¬ 
ticle to scatter. This assumption is valid as At approaches zero or 
numerically where At is chosen to be much smaller than needed 
to avoid multiple collisions. At each time step our numerical code 
calculates the relative velocity and density p using the same func¬ 
tional SPH kernel as in hydro calculations of each DM particle in 
relation to its 32 nearest neighbors. These values are used to cal¬ 
culate the probability that the DM particle interacts with one of its 
neighbors. Similar to Vogelsberger et al. (2012), a scattering event 
may occur at each time step between particles i and j if a uniform 
random variable in the interval (0,1) is less than Vij /2 (the prob¬ 
ability is divided by two in order to account for the fact that each 
set of particles are compared to each other twice). While producing 
similar results, this approach differs slightly from the one adopted 
in [Rocha et al.| ( |2013| ), where collisions are defined between DM 
particle pairs and the probability of interaction P(i|j) is explicitly 
= PG|i)- When a particle collision is detected we isotropically and 
elastically scatter the particles to random angles. 

The interaction cross section a for all SIDM runs in this work 
was set to the relatively large value of 2 g/cm^, close to the up¬ 
per limit for constant cr, in order to maximize its effects compared 
to that of baryon physics. The way in which SIDM reshapes dif¬ 


ferent halos can be readily understood setting N — 1 in eq. 1 and 
V — Umax (where Umax = mux sjGM[< r)/r) to obtain a char¬ 
acteristic maximum timescale on which a given density is stable, 

TSI = -- . (2) 

pUmax O' dm 

After a collision, the change in momentum is usually enough to 
fully eject the particle ( [Kahlhoefer et al.[2014[ l, so dark matter den¬ 
sities must drop over a few rsi. 

One important conclusion to draw from eq.3 is that not every 
SIDM model will necessarily form significant cores in the smallest 
halos. Low halo peak velocities and low central densities (at a fixed 
radius) may result in rsi comparable to the lifetime of a halo and 
therefore in the preservation of cusps. Another important conse¬ 
quence is that in SIDM the central densities of very small satellite 
halos could be differentiated from those of their held counterparts 
by the interaction with the DM halo of a more massive host as the 
orbital velocity is much higher than the internal velocities of the 
satellite. We test these analytical predictions in the next section, 
where we present results of simulations that resolve the internal 
structure of DM halos smaller than most previous SIDM studies. 
Observational evidence of DM cores in small held halos where 
baryonic processes may be inefficient (Vmax < 20-30 km s“^and 
in satellite galaxies (which move at 100-200 km s“^in a ~ 200- 
1000 X Peril DM held) can then provide useful lower limits to the 
SIDM cross section. 


3 CDM AND SIDM SIMULATIONS 

In all simulations we assumed a A dominated cosmology (Do = 
0.26 A = 0.76, (78 = 0.77, n=0.96) and used the “zoom-in” tech¬ 
nique to achieve high resolution in the regions of interest ( [Katz &[ 
[White] 199^ . The gravitational force spline softening length is in 
the range 64-86 pc ( [Power et al.|2003[ l and the smoothing length 
for the gas component is allowed to shrink to 0.1 the force soften¬ 
ing. Simulations start at z = 120 — 100. The combination of sim¬ 
ulations (see Table 1) allow us to compare the effects of SIDM and 
baryon physics on a range covering almost four orders of magni¬ 
tude in halo masses and a variety of environments from the held to 
the dense region within the virial radius of a large held disk galaxy. 
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Figure 1. DM-only Simulations: The DM density measured at 500 
pc as a function of halo maximum velocity and dark matter model. Blue 
crosses: CDM, red crosses: SIDM. Circled symbols correspond to isolated 
halos, non-circled symbols show the satellites of the MW-sized, hl48 halo. 
At small halo masses the central densities of CDM and SIDM halos do 
not differ substantially, while at larger velocities/halo masses SIDM-only 
halos have lower central densities than their CDM-only counterparts. This 
difference can be understood in terms of the dotted lines which show the 
maximum timescale tsi at which collisions are significant in the SIDM case 
(see text for details). This result shows that fixed cross sections commonly 
adopted at the scale of larger systems (0.1-2 cm^/g) would not be sufficient 
to form kpc sized cores in the smallest observable field halos with Yrnax < 
30 km s“^. In the Ymax 30-60 km s“^range SIDM satellites have central 
densities lower by a factor two compared to their field counterparts, due 
the added boost to SIDM collisions coming from satellites orbiting at high 
velocity in the dense DM halo of the host. 


Cosmological simulations of well resolved halos with mass 
< 10^° Mq (corresponding to a halo peak velocity Ymax < 30 
km s“^) are particularly relevant for SIDM, as baryonic effects at 
these scales should be limited ( [Governato et al.|2012[ [Penarrubi^ 
let al.|2012| ). Based on equation (2), the ability of SIDM to lower 
the central density of halos depends on the DM p and the typical 
particle velocity, which is lower in less massive halos. Interestingly, 
previous CDM simulations have suggested that the central density 
of small mass field halos decreases with their halo mass ( | Avila- 1 
|Reese et al.|2005|[Li et al.|2008| ). In a SIDM scenario, extending 
this trend to very small halos could prevent the formation of low 
density central cores by increasing the time scale of DM-DM colli¬ 
sions. 


3.1 DM-only Simulations: Environmental effects and cuspy 
halos at Ymax < 30 km s“^ 

To study the structure of a significant sample of DM halos in a 
range of environments we simulated four different regions with the 
zoomed-in approach ( |Katz & White|1993] ). The first three regions 
(Table 1) are centered on filamentary structures with a density close 
to average measured inside a sphere of 5 Mpc in radius. The largest 
one of them is nicknamed The Forty Thieves’ and includes sev¬ 
eral tens of halos with Ymax < 30 km s“^, equivalent to a mass 
range lO^-lO^^M©, where the minimum mass refers to halos with 



Figure 2. The relative DM content as a eunction oe radius in 
DM-only Simulations: The ratio of the cumulative DM content mea¬ 
sured for a set of small halos over the 0-1.5kpc radial range (from the 40- 
Thieves run). Each CDM halo was matched with its ’twin’ halo in the SIDM 
simulation. The spline kernel force resolution is 64pc. Vertical lines mark 
150 and 500 pc (2.3 and 5 softening lengths). With our adopted cross sec¬ 
tion of 2 cm^/g small CDM and SIDM become progressively difficult to 
differentiate, with DM content at a given radius being within a factor of two 
in the two cosmologies. Detecting such a difference in field dwarfs could 
require a next generation telescope (M.Walker, private communication). A 
SIDM with a variable cross section could show a more significant difference 
compared to CDM. 


at least 50,000 DM particles within the virial radius. To comple¬ 
ment this dataset, a high resolution simulation of a massive halos 
and its system of satellites has been included (hl48). The halos in 
this simulation have mass and spatial resolution similar to their ‘40 
Thieves’ counterparts and mass and spatial resolution better by fac¬ 
tor two compared to recent work |Vogelsberger et al.|p01^ ; |Zavala| 
|et al.| ( [20T^ . In Appendix 1 we verified that our results and the 
Monte Carlo implementation of SIDM collisions are not affected 
by resolution effects. 

In Figure we show the density measured at 500pc (a radius 
at which the rotation curves of many dwarf galaxies have been re¬ 
solved, e.g., |Oh et al.|201l] ) for all resolved halos in our hi48 (x) 
and “the 40 Thieves” (+) DM-only simulations as a function of 
Vmax- We then over plot lines of fixed rsi for 1 Gyr and 5 Gyr. 
SIDM and CDM show identical central DM densities if the typical 
scale for DM interactions is longer than 5 Gyr. At larger peak ve¬ 
locities (and halo masses) with shorter timescales for interaction, 
the SIDM densities decrease compared to their CDM counterparts 
and fall in a valley between the 1 and 5 Gyr lines, matching re¬ 
sults from previous works ( [Rocha et al.|2013[ jVogelsberger et al.] 
|2012| ). In other words, SIDM interactions force the DM density to 
drop until, going to smaller and smaller halo masses, the interac¬ 
tion timescale rises to a significant fraction of the Hubble time and 
SIDM (with a cr = 2 crc?/g) is not effective anymore. Figure]^ 
shows the SIDM vs CDM ratio of the enclosed mass as a function 
of radius for paired halos from the ‘40 Thieves’ run. At radii cur¬ 
rently tested by observations (100-500pc), but entirely neglecting 
baryonic physics, the difference between the two models decreases 
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Figure 3. Radial Profiles. Left:The radial density profile of the DM component for halos/galaxies h516 and h2003. RightiThe velocity dispersion of DM 
component in the same halos/galaxies. h516-CDM: blue. h516-SIDM: red. h2003-CDM: cyan. h2003-SIDM: magenta. Dashed: DM-only runs. Both density 
and dispersion profiles become similar in CDM vs SIDM once baryonic processes are introduced. 



Log halo mass (MG) 

Figure 4. Stellar Mass/Halo Mass relation: The Stellar Mass - 
Halo Mass relation for the simulated galaxies. Dashed lines and solid lines 
show the relation obtained fro m Local Group data jBrook & Di Cintio|2014| 
[Garrison-Kimmel et al.|2014} The relations are extrapolated below ~ 10^'^ 
MQdue to sample incompleteness (dotted lines). Circles show the raw data, 
solid dots show the simulation data correcting for observational and sim¬ 
ulation biases jMunshi et al.|2013) in measuring stellar and halo masses. 
Overall the simulations produce the right amount of stars. The most mas¬ 
sive halo converts about 1% of gas into stars. The rapid drop in SF efficiency 
at halo masses below 10^° MqIs due to the introduction of ’early feedback’ 
(see text). 


rapidly. A factor of two in the DM content within 100-500pc of 
faint field dwarfs, could be revealed by large spectroscopic samples 
of the associated stellar populations (M.Walker, private communi¬ 
cation). We also verified that at larger Vmacc the density profiles of 


the D M-only halos match those recently published in ( [Elbert et al.| 
2014] ), which at Ymax ~ 35 km s“^form significant cores (see also 
next section where the effect of baryons in these larger systems is 
included). This result confirms our simple analytical expectations 
and shows that even with a significant constant SIDM cross sec¬ 
tion, DM cores rapidly become smaller than our resolved scale (~ 
100-200 pc or twice the spline kernel softening) in field halos with 
virial mass < 10^° M 0 and Vmacc < 30 km s“^. Confirming the 
results from [Li et al.j ( |2008| ) and extending them to much smaller 
halo masses and higher resolution we verified that the increased 
Tsi at small halo masses comes not only from a lower aoM but 
also from lower cusp densities, possibly due the later epoch of col¬ 
lapse of the central regions of a halo. We verified the DM density 
within the central 250 pc decreases by a factor of eight in CDM 
halos with mass from 10® M© compared to halos of 10^° M©. If 
we look at the average densities as a function of Vmacc, another dif¬ 
ference emerges between the SIDM field (circled red crosses) and 
satellites (red crosses), with satellite halos showing central densi¬ 
ties lower by about factor of two compared to field halos of similar 
Vmacc - This important environmental difference could be due to the 
satellites forming their central regions earlier, or to the significant 
boost to p and Wmax in eq.3 due to their orbiting 1) inside the dense 
halo of a much more massive host and 2 ) at a much higher speed 
than their internal velocity dispersion. Our simulations were able to 
highlight this difference as our simulations resolve field and satel¬ 
lites halos with significantly lower peak velocities ( 10-20 km s“^) 
than in previous works. 

The main result from this section is that evidence of DM cores 
in real galaxies with Nmax < 30 km s“^ would constrain the SIDM 
cross section to values a »2 cm^/g when the typical DM velocity 
dispersion is low. If significant DM cores are found in these galax¬ 
ies, their existence would give support to models with a variable 
SIDM cross section that is higher (20 cm^/g, see eq.2) at small halo 
masses and then declines rapidly at the scale of groups and galaxy 
clusters as constraints at scale support small cross sections a < 
1 cm^/g. We plan to further explore the relative effect of stronger 
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Run Halo Mass 
ID CDM/SIDM 


V. 
kms 


max 
-1 


Stellar Mass 
CDM/SIDM 


HI mass 
CDM/SIDM 


h516 4.1 xlQio 

h516b 1.1 xlQio 
h2003 1.1 xlQio 


58.4 
33.0 

30.5 


6.2/4.9 X10^ 
9.8/14.7 x 10® 
3.1/3.4 xlO® 


3.3/2.0 xlO®M0 
1.6/0.94 xlO®Mo 
5.0/16.6 x1O8M0 


Table 2. The total halo mass, stellar and HI masses and \max 
(in km S“^) for some representative halos in the CDM runs with SF and 
their SIDM counterparts. z=0 stellar masses were measured within 2.5kpc. 
h516b is the second most massive halo in the h516 run. 




0 


^ Time (Gyrs) 


10 


Figure 6. The slope a ol the DM density proeile in SIDM and 
CDM OVER DIEEERENT MASS RANGES: blue and Cyan: CDM, red and 
magenta: SIDM, solid lines: baryon+DM runs, dashed lines: DM-only runs. 
Top: halo h516. The DM slope evolves rapidly and in a similar way in both 
CDM and SIDM, as dynamical heating is very efficient at this scale. Middle: 
halo h2003. (blue: CDM, red: SIDM). Bottom: a collection of field halos 
with total mass <10® Mq from our simulations. SF efficiency is too low 
and the SIDM interaction rate is too low to create cores. 


Figure 5. Eeeects oe Baryons Physics: The central density of CDM- 
only halos (blue) vs their SIDM (red) counterparts run with baryon physics 
and bursty feedback. Lines connect the DM-only runs of h516 and h2003 
to their counterparts with baryons. For our choice of the SIDM interaction 
cross section and for halos with Nmax > 50 km s“^ (corresponding to stel¬ 
lar masses larger than 10® Mq) the central mass of CDM and SIDM halos 
is similar, but significantly lower than the predictions of CDM-only runs. 
Due to low DM interaction rates and lack of bursty outflows, at lower halo 
masses both cosmologies give similar predictions: cuspy central DM pro¬ 
files. Larger SIDM cross sections would be able to differentiate the various 
models in halos with Ymax < 20km s“^as DM cores would then form in 
galaxies were baryon physics do not play a major role. 

SIDM interactions on satellites compared to field halos in future 
work. 

3.2 Simulations with SF: SIDM similar to CDM when 
Baryon Physics are relevant. 

We focused our hydrodynamical simulations on the largest galax¬ 
ies formed in the filamentary regions ‘h516’ and ‘h2003’. These are 
two well studied fields: h2003 is the same halo simulated in IGov- 
|emato et ari ( |2015| ) while h516 is the main halo of |Governato et al. 
( 2010 1 and of the “7 dwarfs” galaxies sample studied in |Shen et~^ 
(20141 and |Madau et al.| ( [2014| ). The SF parameters in our study 
were identical for all CDM and SIDM simulations. They also cor¬ 
respond to the fiducial “g5” runs in |Govemato et al.| ( |2015| ), where 
we explored the effects of different feedback and SF recipes in the 
context of comparing the formation of dwarfs in CDM vs Warm 
scenarios. Here we emphasize that our SF implementation creates 
repeated starburts and gas outflows with significant loading factors 
(gas mass ejected from the center divided by star formation rate). 


Multi-wavelength evidence for outflows, analysis of the stellar pop¬ 
ulations in the SDSS dwarfs and realistic CMDs ( |Governato et al.| 
|2015| ) give strong support to our implementation of SF. This ap¬ 
proach differs from the SIDM study of |Vogelsberger et aL] ( |2014| ) 
where less bursty feedback still removes gas from galaxy centers, 
but does not create substantial DM cores. The baryonic content of 
the galaxies in our study are summarized in Table 2. 

Figure (left panel) compares the density profiles of the DM 
component of galaxies h516 and h2003. In DM-only simulations 
halos with masses > 10^° Mq have cuspy profiles (halo h2003: 
cyan dashed, halo h516: blue dotted). Once baryon physics and 
outflows are introduced, flatter DM profiles are created in both 
SIDM and CDM cosmologies. The blue dashed (CDM-only) vs 
blue (CDM-i-SF) lines and the red (SIDM-i-SF) show results for 
h516. The cyan dashed (CDM-only) vs cyan solid (CDM-i-SF) lines 
and the magenta (SIDM-hSF) lines show results for h2003, the 
smaller halo. In both CDM and SIDM models the central cuspy pro¬ 
files have been significantly flattened inside 1 kpc. h516, the most 
massive halo studied with the inclusion of SF, is the one where 
the central DM density decreases the most. This result is consis¬ 
tent with previous findings ( |Governato et al.|2012||Di Cintio et al!] 
|2014| ), showing that the efficiency of core formation peaks in halos 
with Ymax ~ 50 km s“^ 

Figurej^ (right panel) compares instead the velocity dispersion 
profiles of the DM component in halos h516 and h2003. As for the 
density profiles, DM-only runs of different cosmologies have dif¬ 
ferent local DM velocity distribution profiles: CDM halos show a 
decreasing dispersion closer the halo center while the SIDM ha¬ 
los show a rather flat profile. This difference, a result of the energy 
transfer to the center of the halo due to collisional processes, had 
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Figure 7. The projected color density map op the baryon distribution in galaxy h516. Top: CDM. Bottom: SIDM. Left panels: projected 
gas density at z=0 seen edge-on and face-on. Right panels: Stellar projected density at z=0. The stellar disks are relatively thicker than in earlier simulations 
( jGovernato et al.polo) , an effect of the introduction of ‘early feedback’ from young stars. We verified that in both models the stellar distribution is exponential, 
with a similar disc scale length and lack of a central dense spheroid. 


been considered a strong signature of SIDM. However this differ¬ 
ence is erased by the introduction of bursty feedback, that creates 
significant DM cores in the CDM halos. ‘Dynamical heating’ also 
causes the velocity dispersion profiles of all halos to flatten out. 

Figure shows the stellar mass/halo mass ratio of the halos 
we simulated with the inclusion of baryon physics, showing that 
they follow the stellar mass/halo mass relation inferred using local 
data and the abundance matching technique ( [Brook & Di Cintio| 
|2014[ [Garrison-Kimmel et al.||2014j ). This is particularly relevant 
as producing the correct amount stars is necessary to estimate the 
minimum mass at which baryonic processes can originate cores. 
After this paper was submitted ( jOnorbe et al.|2015] ) published re¬ 
sults from a simulated CDM halo of similar total mass, which also 
formed abut 10® MqIu stars, leading to a cored DM profile. 

Figure]^ shows the DM central densities as a function of halo 
peak velocity once SF is included. The comparison with Figure[^is 
striking: where at Vmacc > 30 km s“^CDM was clearly differenti¬ 
ated from SIDM, now the central density of CDM and SIDM halos 
is almost identical over the whole 10 to 50 km s“^range. This result 
clearly illustrates how predictions form DM-only runs can be com¬ 
pletely superseded by the addition of the complex baryon - DM 
interactions. 

In Figure we investigate how the slope of the DM profile 
(measured at 500 pc) evolves with time in a sample of field halos 
with a range of masses. As the mechanisms of core formation dif¬ 
fer (potential fluctuations linked to SF vs SIDM collisions) the time 
evolution of the DM slope may also be significantly different. We 
have found that the evolution of the systems can be roughly divided 
in three stellar mass ranges. Present day systems with stellar masses 
larger than 10® M© (Ymax > 50km s“^as in classic field dwarfs) 
had cored DM profiles since redshift > 4, with no significant differ¬ 


ence between SIDM and CDM (top panel). In the SIDM-only run, 
the DM core forms more slowly than when the effect of baryons are 
included. The effect of baryon feedback on the DM profile is clearly 
detectable as soon as about 10® M© of stars have been created (see 
also jGovernato et al.|2015[[Fenarrubia et al.|2012j l. In galaxies with 
present day stellar masses below 10^“® M©, the effect of ‘dynam¬ 
ical heating’ induced by feedback is progressively reduced ( jGover-j 
jnato et al.|2012|[20T5l|Di Cintio et al.|2014| l. At the radial scale of 
500 pc the CDM and SIDM halos of h2003 show DM profiles that 
start cuspy and then slowly develop a flatter profile. In both SIDM 
and CDM the process is more gradual compared to more massive 
systems, due to a number of factors (star formation rates are lower 
and Tsi is longer). However, the inclusion of baryonic processes 
again makes the galaxy evolve similarly in CDM vs SIDM. By red- 
shift z<l the difference in DM slopes is not significant enough to 
discriminate between SIDM and CDM. Halo h2003 was also re¬ 
cently run in a WDM cosmology ( jGovemato et al.|2015j l. Similar 
to the CDM case the WDM cusp turned into a core over time, due to 
bursty feedback. By z=0 the central DM distribution of halo h2003 
is similar in CDM, WDM ( jGovernato et al. 12015) 1 and SIDM (this 
work). 

At even smaller masses, the evolution of the DM slope a with 
the inclusion of baryonic processes (continuous lines) confirms the 
results from the DM-only runs (dashed lines). SF at such small 
scales is strongly inhibited by the cosmic UV field. Similarly, at 
such small masses the SIDM interaction timescale becomes long 
compared to the Hubble time and the survival time of the halo. 
Hence the central DM slope measured at 500pc is not significantly 
different in SIDM vs CDM halos, as predicted by the empirical cal¬ 
culations in the previous section. The failure of SIDM to form sub¬ 
stantial DM cores in very faint dwarfs with host halos with mass 
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Figure 8. The SFR as a function of cosmic time. Left panel: halo h516 in CDM and SIDM cosmologies. Both galaxies have extended, gas rich disks 
and negligible bulges. Right panel: the burstiness of models h516 and h2003 in CDM and SIDM, measured at different timescales (see text for a definition). 
SF does not differ significantly in CDM vs SIDM, and appears to be mostly regulated by feedback and accretion rates. 


< 10^° M© had not been explored in previous numerical works, 
where high resolution simulations mostly focused on halos more 
massive than 10^° M©. 

Overall, the inclusion of SF and baryonic processes erase the 
differences seen in the DM-only runs. At the scale of dwarf galax¬ 
ies with Vpeak ~ 30-50 km s~^. energy transfer to the DM makes 
the DM distribution at the center of galaxies flatter than an NFW 
profile and almost identical in CDM vv SIDM . As a result cur¬ 
rent constraints on the SIDM cross section from dwarf sized field 
come out considerably weakened and will have to be reconsidered. 
Moreover, we predict (see Figure that a SIDM cross section > 
10 cm^/g would necessary to create DM cores in smaller halos with 
Vmacc < 10 km s“^. As the shapes of gal axy clusters requir e a 
SIDM cross section smaller than 1 cm^/g |Peter et ^|2013| ). If 
cores exist in dwarf galaxies, these results argue in favor of a vari¬ 
able SIDM cross section (see also |Vogelsberger et al.| ( |2012] )). 


3.3 Star Formation Histories and Assembly of the baryonic 
component 

In this section we analyze the z — ^ baryonic distribution of the 
galaxies that form stars in our simulations, while focusing on the 
two main halos: h2003 and h516. We verified that in both CDM and 
SIDM cosmologies the stellar mass of galaxy h516 follows an ex¬ 
ponential profile with no spheroidal component. The projected gas 
and stellar densities of galaxy h516 are shown in Figure show¬ 
ing a similar distribution in CDM vs SIDM. Both stellar discs are 
thick and dynamically hot (Figure and more extend ed than the 
CDM version of h516 described iii |Christensen et ^ ( 2014^ , an 
effect of the introduction of early feedback ( |Trujillo-Gomez et'^I] 
|2013[[Ruskar et al.|2014| t. Confirming the visual impression from 
Figure [7] we verified that the radial distribution of stars at z=0 does 
not change significantly in the CDM galaxies vs their SIDM coun¬ 
terparts, both following a trend for larger systems to have more 


extended stellar systems and roughly consistent with observational 
data. The half light radius for h516b (Table 1) is for example, 2.0 
kpc in SIDM, vs 1.8 kpc in CDM. This is a different outcome com¬ 
pared to the results of a less bursty SF scenario as the one studied 
in |Vogelsberger et al.] < |2014| ), which followed galaxies in a simi¬ 
lar mass range to ours. In their simulations the SIDM distribution 
is less concentrated compared to the CDM one, leading “to dwarf 
galaxies with larger stellar cores and smaller stellar central densi¬ 
ties.” In our simulations the dynamical coupling of baryons to DM, 
driven by fast and repeated outfiows, shapes the stellar distribution 
and leads to a similar baryon distribution and content. 

Furthermore, with the adopted SF prescription, the star forma¬ 
tion histories (SFH) of the galaxies formed in the CDM model look 
very similar to the SIDM counterparts. This is most likely due to 
two main reasons: 1) the assembly rate of both DM and baryons, 
being driven by the large scale structure remains unchanged in 
SIDM and 2) there is no difference in the central DM distribution 
between the two models as it is responds quickly to the effects of 
feedback. As an example. Figure shows the SFH of halo h516 in 
the CDM and SIDM models, which show a very similar time evolu¬ 
tion, peaking 5 Gyrs after the Big Bang. Overall, SIDM and CDM 
halos have almost identical DM masses by z=0 and the SF effi¬ 
ciency is similar, within 20% in the two models. This result differs 
from the CDM vs WDM comparison in |Governato et al.| ( |2015| ), 
where the WDM version of h2003 formed only half the stars due 
to the delayed assembly of the halo (see also |Colm et ^|2015| l. 
However, changes on the timescale and intensity of the individual 
SF events could be driven by subtle effects and still differentiate 
between CDM and SIDM, for example a shallower density profile 
could drive more gas instabilities that lead to radial inflows and in¬ 
crease the burstiness of SF, or, on the other hand, lower DM densi¬ 
ties could slow down the collapse of gas and reduce SF in the SIDM 
scenario. We quantified the burstiness of SF in SIDM vs CDM by 
measuring the dispersion in star formation rates (SFR) measured 
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Figure 9. The rotation curves Vc of halos h516 and h2003 in CDM (blue and cyan) and SIDM (red and magenta) compared with the observational 
constraints from a sample of Local Group field dwarfs (Wolf et al.|2010|[McConnachie|2012||Weisz et al.|2014) . The introduction of feedback processes 
lowers Vc out to at least 2 kpc, where mass decomposition based on HI kinematics are crucial. (Vc is defined as = (GM/r)^/^). Softening is 64pc for h2003 
and 86pc for h516. 


in AT SFR intervals, where AT sfr was varied from 10® to 10^ 
years (for a similar approach but different definition of burstiness, 
see |Hopkins et al.|2014| l. Measuring AT sfr over hundreds of ran¬ 
domly sampled intervals, bursty SF at a given time scale shows as a 
large dispersion if the SF rate at a time T+AT^ fr differs substan¬ 
tially from the SFR at a previous time T. The right panel of Figure 
m shows that the burstiness measured over different timescales was 
similar for both halo h516 and h2003 in both models. Both h516 
and h2003 show a SF rate dominated by small time scale fluctua¬ 
tions, with a similar dependence on time scale and halo mass, with 
the less massive halos having a burstier SF at all time scales (a re¬ 
sult seen also in simulations with different feedback models, e.g., 
[Hopkins et al.||2014) . Observations of nearby stellar populations 
and comparisons with simul ations ([Tolstoy et al.|200^|Kauffmann| 


[2014[[Governato et al.|2015[|Weisz et al.|2014 l strongly support a 

bursty build up of the stellar content of dwarf galaxies. We verified 
that our set of simulations qualitatively reproduce the fraction of 
stars formed in bursts as estimated in [Kauffmann[ (2014] )): 20-50% 
with a larger fraction formed in the less massive system. Bursts 
were defined as periods with SF > 2 —4x < SFH > and duration 
of 10^-10® years. A burstier implementation of our SF model (as 
one that neglected HI self shielding) would be ruled out by exces¬ 
sive structure in the stellar CMD, as discussed in [Governato et al.| 
( |20T5l ). 


Overall, Within our small sample of galaxies, we do not find 
evidence of large differences in the baryonic content and distribu¬ 
tion in galaxies formed in a CDM SIDM model 


Finally, Figure show the rotation curves Vc (defined as 
^JGM{< r)/r) of our galaxies h516 and h2003 compared to es¬ 
timates from a sample of local dwarfs (empty squares, from [Wolf[ 
[et al.[2010[[Weisz et al.[2014] ). The left panel shows the well known 
result that halos formed in CDM-only simulations have higher cen¬ 
tral mass densities that most observed dwarfs (although significant 


uncertainties remain on the observational estimates). At the same 
time, SIDM-only runs with a range of cross sections (Elbert et"^ 
[2014[ ) show a better agreement with observational estimates as the 
central DM densities are lower. However, the right panel of Figure 
1^ shows that once baryon physics are introduced, the central mass 
distribution (DM+baryons) of CDM dwarfs matches the observa¬ 
tional data, equally well, without requiring SIDM. Recent works 
have compared dwarf galaxies in the Local Group (Brook & Di 


[Cintio|20T4 Sawala et al.|2015[ ) and MW satellites [Zolotov et al. 
(2012[ |; [Arraki et al.[ (2014[ ) with simulations. They have shown 
how the introduction of feedback (and tidal processes for satellites) 
generates a realistic M^tar - ^haio relation, as galaxies at a fixed 
Mstar are hosted by more massive halos than when a universal 
NFW profile is assumed in the analysis This result shows how 
the perceived discrepancy between the observed number density of 
field galaxies as a function of Wmax and that of the underlying halo 
population is naturally accounted for in a CDM context where feed¬ 
back is introduced and the central total density is lowered. In good 
agreement with results from our simulations, the revised abundance 
matching from [Brook & Di Cintio[ ( |20T^ predict that field galax¬ 
ies are hosted in halos with mass >5 x 10^ Mq, with most ha¬ 
los below that threshold being devoid of stars. Our results provide 
direct evidence that SIDM-only simulations mimic the effects of 
feedback, but also that when bursty feedback is introduced, SIDM 
galaxies are essentially indistinguishable from their CDM coun¬ 
terparts. Both SIDM and CDM models create galaxies with mass 
distributions and observable properties in broad agreement with the 
observed ones. 


^ in a halo model with a central mass density lower than NFW the Vcirc 
measured at ~ 1 kpc translates to the same Ypeak and hence halo mass of 
an NFW halo. As a consequence, an observed Wcirc translates to higher 
halo host mass than when a more concentrated NFW halo is assumed. 
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4 CONCLUSIONS 

We Studied a set of cosmological simulations of field dwarf galax¬ 
ies in CDM and SIDM. The SIDM simulations include for the first 
time a description of bursty SF and feedback that creates potential 
fiuctuations able to lower the central DM density of halos where 
at least 10® M© of stars have formed. The force resolution of our 
simulations (64-86 pc with a spline kernel softening) allows us to 
resolve with many thousands of particles the central regions of ha¬ 
los (and millions to the virial radius) down to halo masses where 
the impact of star formation will not play a major role and where 
the effect of SIDM alone dominates. A relatively large value (2 
cm^/g for the SIDM cross section was chosen to maximize its dif¬ 
ference from the CDM model. We have also run a set of DM-only 
simulations that allowed us to study the evolution of a uniform res¬ 
olution set of DM halos down to masses of only few times 10® 
M© (equivalent to Vma® < 20 km s“^) both in the field and in 
the dense environment of a more massive host. Our results can be 
summarized by three main findings: 

• Once SF and resulting feedback is introduced, the central DM 
mass distribution and velocity dispersion becomes similar for CDM 
and SIDM galaxies with stellar masses in the 10®-10® M© range 
and circular velocity 25-50 km/sec (Fig.3). At the scale of 0.5-2 
kpc the total matter content is in good agreement with observa¬ 
tional estimates of Local Group dwarfs. This result differs starkly 
from the predictions of DM-only simulations, or simulations where 
energy feedback does not lead to core formation ( |Boylan-Kolchin| 
|et al.|20lT||Garrison-Kimmel et al.|20f4] ). 

• Analytical expectations and simulations show that with a fixed 
SIDM cross section of 2 cm^/g the DM central density and internal 
velocities of halos with mass below 10^ M© are too low to have 
a significant number of DM-DM interactions. Such SIDM halos 
remain cuspy when observed at a scale of 500pc, and show an en¬ 
closed DM mass content lower by about only a factor 2 compared 
to their CDM counterparts (Fig.l & 2). This result poses an inter¬ 
esting lower limit to the SIDM cross section in halos associated 
with galaxies with stellar masses below 10® M©. If kpc sized DM 
cores are found in these galaxies, their existence would give support 
to models with a variable SIDM cross section that is high ( ~ 10- 
20 cm^/g, see eq.2) at small halo masses and then declines rapidly 
at the scale of galaxies and clusters. Large cross sections at dwarf 
scales are also suggested by DM-only studies ( [Elbert et al.|2014| |. 
Such large cross sections may have a detectable effect on the SFHs 
of dwarf sized systems, but could also make galaxy satellites easier 
to disrupt. Our simulations also highlight the relative environmental 
effects of the enhanced interaction rate at the center of very small 
satellites compared to their field counterparts, due to rapid orbital 
velocities in the dense halos of a massive host. This relative dif¬ 
ference offers the potential of using faint galaxies to constrain the 
cross section of SIDM as a function of the DM velocity. 

• Once regulated by feedback, the SFHs (Fig.7), stellar and gas 
content and spatial distribution do not differ substantially in CDM 
vs SIDM model, both forming gas rich galaxies with bulgeless 
disks which resemble real ones. 

These results lead to potentially important implications for 
SIDM models. As baryons lead to a drastic change in the mass 
distribution of central regions of typical dwarf galaxies with Vc 
~ 30-60 km/sec, many of the current bounds on the SIDM cross 
section at different mass scales will need to be re-examined. After 
the introduction of realistic descriptions of SF and feedback. Both 
CDM and SIDM form dwarf galaxies with DM distribution and 


observable properties in broad agreement with observational data. 
As a result, while certainly not excluded, SIDM is not necessary to 
solve the problem of the excess of DM at the center of simulated 
dwarf galaxies with Wmax > 30 km s“^. Furthermore, because 
the effect of SIDM becomes negligible at very small galaxy/halo 
masses even with a relatively large constant velocity cross section 
(2 cm^/g), SIDM with a velocity dependent cross section will be¬ 
come necessary if even very faint dwarf galaxies are observed to 
have dark matter cores. This work highlights how simulations that 
a) include realistic SF and feedback processes and b) study a large 
sample of faint dwarfs are necessary to take advantage of astro- 
physical constraints on DM models. 
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Figure 10. The density in scaled units of R/Rs for Hernquist halos of 
10^° (black), 10^ (green), and 10® MQ(blue) which have scale radii of 
1.099, .37, and .12 Kpc respectively based on the mass concentration ratio 
relation from Neto et al 2007. Each halo has been evolved for about 25 
scaled dynamical times which corresponds to 91, 44, and 20 million years 
respectively. The halos demonstrate self similar evolution in these scaled 
distance and time units. The density p has been scaled by the central density. 
The vertical dashed lines indicate the softening length of 7.5 pc that is scaled 
for each halo (the dashed line for the 10 ^° M 0 halo is off the plot further to 
the left). The dashed grey dashed line is the Hernquist density profile. 
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Figure 11. Snapshots of a Hernquist profile SIDM simulation showing the 
number of interactions that occur per particle per giga-year for the same 
halos as above. The cross section is 10 cm^ g“^. In dashed grey is the an¬ 
alytical prediction for each halo respectively. The simulations were run for 
a period much shorter than the dynamical time of the halo. Each Hernquist 
profile has a different scaled structure as described in the previous figure. 
Each halo was described with about 2 million particles. 


an integral of the velocity weighted cross section in the volume V. 

r = J ^;^{admVr)ix)dV (3) 

Where p is the local DM density, is the mass of the dark matter 
particles in the simulation, Udm is the dark matter cross section 
per gram, and {admVr) (x) is the local thermal average of the cross 
section weighted by the particle’s relative velocity as a function of 
the position x. The local thermal average of the cross section is 
calculated by taking the first moment of admVr over the combined 
relative velocity distribution function f{vr). 


5 APPENDIX: THE SIDM IMPLEMENTATION 

We verified our implementation of SIDM by comparing analytical 
predictions of the number of dark matter interactions to simula¬ 
tions of dark matter halos with a Hernquist density profile ( |Hern-| 
|quist|1990] l (Fig.[^. We then demonstrate the formation of a flat 
core and the signature flat SIDM velocity dispersion profile, start¬ 
ing from an NFW profile with a cuspy density profile and a rais¬ 
ing velocity dispersion ( [Navarro et al.|I996] l. Our implementation 
code closely replicates published results of the predicted number 
of SIDM collisions as a function of halo mass and local density 
and the evolution of the central density profile ([Vogelsberger et al.j 
[20T2l[Rocha et al.|20I3[ l (Fig[^and Fig[TTJ. 

The number of dark matter interactions per unit time is an 
important cosmological value for SIDM theories and a useful di¬ 
agnostic for numerical implementations of SIDM interactions The 
probability of an interaction between any two particles is given in 
eq.4. The total number of interactions, F, that occur per unit time is 


(aUr)(x) = / (admVr){x)f{Vr)dVr (4) 

Jo 

In the most general case the distribution function can be a function 
of position (as for a dark matter halo) and the cross section can 
be a function of relative velocity (for velocity dependent SIDM). 
In the simplest case for single value initial velocity of vq and uni¬ 
form density, the distribution function f{vr) becomes a delta func¬ 
tion at that given velocity where f{vr) = and the inter¬ 

action rate simplifies to F = ^/2l2NpadmV where N is the to¬ 
tal number of dark matter particles in the simulation. A standard 
Maxwell-Boltzmann distribution has a constant velocity dispersion 
so that after averaging over all angles of collision we conclude 
(Vr) = V2{V) SO that if a is the standard deviation of the velocity 
vector in the distribution then the thermal average cross section is 


) (x) — ^^/^CFdm\ (5) 

V TT 

Thus in a constant density region with a Maxwell-Boltzmann ve- 
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locity distribution the number of interactions per unit time is 

V2Npadm tia? 

r =-^(6) 

We would like to have an analytic calculation for the number of 
interactions that would occur in a given dark matter halo, how¬ 
ever in general the integrals of dark matter halos are not well be¬ 
haved. Further, the dark matter interactions modify the density and 
velocity distribution of the halo by conducting energy (as a mea¬ 
sure of velocity dispersion) from the outer regions of halo into 
the core. We approximate that after 25 scaled dynamical times 

— admTn^Rs \J 64Gpq (where po = M j {27rRl) is the central 
density) a core-collapse phase may start ( |Koda & Shapiro|2011| |. 
One well behaved halo profile is the Hernquist profile which has a 
dark matter density profile similar to an NFW profile. In figure 
we show scaled Hernquist halos with masses of 10^°, 10^, and 10^ 
Mq that have each been evolved for 25 scaled dynamical times. 
Under the assumption of an isotropic velocity distribution, the ve¬ 
locity dispersion ay of the Hernquist profile is known analytically 
from the literature and the average of {admVr){'^) is then 

The total number of interactions per unit time is then given 
by equation 4 in the Appendix. However because the initial ve¬ 
locity dispersion is that of a Hernquist profile in equilibrium, this 
analytic prediction is an approximation which declines in accuracy 
with halo evolution in time. In figure we show the analytic pre¬ 
diction of the expected number of DM interactions compared to a 
simulation of the halos that has been run for a brief period ( 10 ^ 
years), much smaller than the dynamical time of the halos and be¬ 
fore the density profile evolves significantly. 

5.1 Convergence of SIDM profiles in DM-only simulations 

In Figurewe show a resolution test using the most massive halo 
of simulation ‘h516’ The halo (of final mass 4 x 10^° M 0 was sim¬ 
ulated lowering the mass resolution by a factor of 64 in mass and 
4 in force resolution. The ratio of the spherically averaged local 
density as a function of radius show that the SIDM simulation con¬ 
verges at about 2 softening lengths, similar to the CDM run (see 
also |Power et al.| ( [2()03| )) The ‘central’ density decreases by a factor 
of ~ 6 between the CDM and SIDM case. These results agree with 
previous works. To verify the convergence of the SIDM density 
profiles at the previously poorly explored regime of M^;^^ < 10 ^ 
Mqwc also simulated the ‘40 Thieves’ volume first with particle 
mass 8000 M 0 and then again 2400 M 0 , and a force resolution of 
65pc in both cases. Combining these runs with the hi48 volume 
we are able to cover four orders of magnitude in halo mass, from 
small halos with peak velocities smaller than 10 km s“^to massive 
galaxies with a large system of satellites. We have simulated hl48 
at lower resolution with hydrodynamics and SF, finding that it hosts 
a large disc galaxy. 

This paper has been typeset from a T^X/ file prepared by the 
author. 
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Figure 12. The local density ratio as a function of radius between the high 
resolution SIDM version of halo h516 (psidm) and its lower resolution 
counterparts (pisse) and (pres) Results in the low res versions converge at 
~ two softening lengths, as typical of CDM-only simulations. 
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